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Abstract 

Using perturbative renormalization group we investigate the influ- 
ence of random velocity field on the critical behavior of directed bond 
percolation process near its second-order phase transition between ab- 
sorbing and active phase. Antonov-Kraichnan model with finite cor- 
relation time is used for description of advecting velocity field. The 
field-theoretic renormalization group approach is applied for getting 
information about asymptotic large scale behavior of the model un- 
der consideration. The model is analyzed near its critical dimension 
through three-parameter expansion in e, 5, rj, where e is the deviation 
from the Kolmogorov scaling, 6 is the deviation from the critical space 
dimension dc and rj is the deviation from the parabolic dispersion law 
for the velocity correlator. Fixed points with corresponding regions 
of stability are determined to the leading order in the perturbation 
scheme. 



1 Introduction 

Percolation processes are the famous models for description of random struc- 
tures [H [2] . Part of these processes are known as directed bond percolation 
(DP). The general feature of DP is that the agent (particle) can propagate 
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from one site to another site in the allowed passage direction. Direction of 
agents spreading is designated by preferred direction of space. The DP can 
also serve for explaining hadron interaction at very high energy (Reggeon field 
theory) [3], various models of spreading disease [H [2], stochastic reaction- 
diffusion processes on a lattice |1] or as in original formulation [5] wetting of 
porous material or exploring path in labyrinth. The upper critical dimension 
for this problem was estimated to be ci = 4 in contrast to the value d = Q for 
the isotropic dynamical case [H El [7] . 

One of the most important property of DP is an exhibition of non- 
equilibrium second-order phase transition between the absorbing and the 
active phase. The absorbing (inactive) phase corresponds to the case, when 
medium does not contain agents (sick individuals) and the active phase rep- 
resents behavior of system, when the number of agents fluctuates around 
constant value. 

A lot of effort was put into the investigation of various effects, e.g. in 
papers [8] long-range interactions by the means of Levy-flight jumps were 
studied both in time and space variables, introduction of immunization was 
examined in [S] , effect of surfaces was studied in (TU] etc. One can also easily 
imagine that spreading of disease can be rapidly enhanced by some external 
atmospheric current or by flying insects. In both cases additional drift can 
be modelled by the means of random velocity field [11] with prescribed sta- 
tistical properties. In this paper the influence of advective field described by 
rapid-change Kraichnan model was studied, which is characterised by white- 
in-time nature of velocity correlator. Generalizing this approach it is possible 
to study e.g. effect of compressibility [T^ or using stochastic Navier-Stokes 
equations effects of "real" turbulent field [13]. In this work we investigate 
the influence of finite correlated velocity field (for introduction see [2]) and 
determine how it can change the critical behavior of percolation process. By 
the means of renormalization group approach we determine possible fixed 
points with corresponding regions of stability. We show that the model ex- 
hibit 12 possible large-scale regimes and the influence of velocity fluctuations 
is illustrated on quantities as the number of particles, the survival probability 
of an active cluster and the radius of gyration of the active particles. 

The paper is organized as follows. In the section 2 we give the brief 
description of the model. In the section 3 multiplicative renormalizability of 
the model is demonstrated and the RG functions (/3 functions and anomalous 
dimensions) are calculated to the one-loop approximation, in the section 4 
possible large-scale regimes are listed and their physical meaning is given in 



2 



the section 5. 



2 Field-Theoretic Formulation of the Model 

Two different approaches for field-theoretic formulation of DP are possible. 
The first approach is more rigorous and is based on the use of master equa- 
tion, that can be rewritten employing Doi formalism [15] into the form of 
time-dependent Schrodinger equation with non-hermitian Hamiltonian. Af- 
ter the continuum limit is performed the effective action is obtained, which 
is amenable to the usual field-theoretic methods. Action for the pure DP 
problem [T] can be written in the following form 

^1 = + DoV'-Doro)iJ + ^[{^^Y^-^^i:% (1) 

where all required integrations over space-time variables are implied. Here 
ip = ip(t,x.) is the coarse-grained density of infected individuals (agents), "0^ 
is the response function, Dq is the diffusion constant, Aq is positive coupling 
constant and Tq is deviation from the threshold value of the injected prob- 
ability. One can assume tq ~ — P, where pc is a critical probability for 
observing percolation (analogous to the deviation from critical temperature 
for equilibrium models). It can be readily shown [1] that the second approach 
based on the use of Langevin equation leads after appropriate rescaling to 
the same coarse-grained action functional ([T]). The both of them lead to the 
same predictions for the universal quantities. 

The agents can be considered as passive scalar quantity [H] that is ad- 
vected by the velocity field with no back influence on the velocity field itself 
with nontrivial interactions given by the cubic terms in ([1]). The inclusion 
of the velocity field v(t, x) corresponds to the replacement partial derivation 
with respect to time by convective derivation = 9t + (v ■ V) in equation 
(II]). In contrast to the work of Antonov et al. [11] we assume that the veloc- 
ity field is a random Gaussian variable with zero mean and correlator in the 
following form [T3] 

(..(t,x).,(t',x')) = 1^1 ^i^?/^.(c.,fc)e-(-*')+^'^(— (2) 

where Pij{k) = 6ij — kikj/k^ is transverse projection operator and Dy has the 
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following form [13] 



^''^"'^) = a;2 + nfo/^o^(P-.)^' 

where 5^10 is the coupling constant and e, r/ play the role of small expansion 
parameters. In this paper e should be understood as deviation from Kol- 
mogorov scaling [16] and 5 is the deviation from the space dimension four 
via relation = 4 — 25. The exponent 77 are related to the frequency uj oc k'^~'^ 
and ?7 = 4/3 corresponds to the Kolmogorov frequency. The averaging pro- 
cedure with respect to the velocity fluctuations may be performed with the 
use of Gaussian-like action functional S2 = —\'vD^^\ . 

The considered model ([2]), ([3]) contains two important limits. The first 
of them is known as rapid-change model [H], that is defined by uiq — )■ 00 
a g'^Q = qio/uIq = const., Dy{u!,k) — )■ g[QDok~'^~'^^~^^ . Hence we see, that 
velocity correlator is decorrelated (white noise) in time variable. The second 
limit is called 'frozen' velocity field [H] and is obtained as Uiq — >■ and 
9o = 9o/uio = const., D^(a;,k) — )■ g'lQDlk~'^~^^~^''7i6{u). In this case the 
velocity correlator is independent of time in the t-representation. 

The full problem is equivalent to the sum of functionals corresponding to 
the DP and velocity field. 

^0 = ^P^[-^t-{^r.V) + DoV'-DoTo]i^+^[{^^f^-^^^']-l^rD;\. (4) 

In perturbation theory it can be easily shown [11], that expansion parameter 
is rather Aq than Aq, which is a consequence of the so-called rapidity reversal 
symmetry (in the language of Reggeon field theory) 

iPit, x) ^ -4j\-t, x), 4j\t, x) ^ -tPi-t, x). (5) 

Therefore the introduction of the following new charge (?2o = -^o appropri- 
ate. 



3 Canonical Dimensions and UV divergences 

Theoretical analysis of UV divergences is based on the standard power count- 
ing [TTIIIH]. Dynamic models involve two independent scales: the time (fre- 
quency) and length (momentum) scale. Therefore the canonical dimension 
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of any quantity F is determined by two dimensions: the frequency dimension 
dp and the momentum dimension dp. These dimensions are found out from 
the usual normahzation condition 

c = -^^r = i, dt = -di = i, dt = -dt = o (6) 

and from the requirement that the action (j4]) should be dimensionless (with 
respect to the momentum and frequency dimensions separately). The total 
canonical dimension is given by the relation dp = dp + 2dp (in order to have 
parabolic dispersion law dt oc for the free theory). It plays the same 
role in the theory of renormalization as length (momentum) dimension for 
static theory [17]. The canonical dimension of the model dl]) are given in 
the Table [H The model is logarithmic at space dimension d = 4 (or 6 = 0) 
and for e = t] = 0. The UV divergences in the minimal subtraction (MS) 
scheme realize themselves as poles in e, 6, rj or their linear combination. The 

Table 1: Canonical dimensions of fields and bare parameters of action 
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Ao 


9w 


920 




U20 
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d/2 


d/2 


-1 


-2 


2 


5 


2e + ri 


26 


V 





d^F 








1 


1 




















dp 


d/2 


d/2 


1 





2 


6 


2e + r] 


26 


V 






total canonical dimension drj^ of an arbitrary 1-irreducible Green function 
Tn = ($...$) is given by the following relation [I71[T8]: dpiv = d+2-iV$d$, 
where A'',^ = {A^„, A^^, A^^t} are the number of fields entering into the Green 
function F^v, including summation over all types of fields $. In logarithmic 
theory superficial UV divergences can affect only those functions F^v, for 
which dpjv is a non-negative integer. Green functions must contain the both 
fields ip and ip'^ because function {tfj . . .ip) or {ip'^ • • • V^^) includes closed loop 
of retarded propagator [HI [17]. Therefore superficial UV divergences can be 
present only in the functions: 

{tp'^tp) with the counterterms ifj^ dtijj , d"^ ijj , ip'^ ip , 
{^"^^ip) with the counterterm V'^V^^; 
{ijj'^ ijj^ ip) with the counterterm ('^p^y^p, 

{ip'^ipv) with the counterterm ip\vd)ip, 
[ip'^ipvv) with the counterterm ip'^ipvv. 
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The terms and ip^ip'^ can be renormalized by the same renormahzation 

constant as a consequence of symmetry and counterterm produced by the 
Green function (tp'^ipv) reduces to the form ■ijj'^{vd)ip = —■ilj{vd)ip^ owing to 
the transversahty of velocity field. In Kraichnan (rapid change) model the 
counterterm ip'^ipvv is absent and renormahzation constants for terms ip^dtil' 
and il)'^{vd)il) are the same (convective derivation Vt = 9t + (v • V) conserves 
its form during renormahzation) due to Galilean invariance of the model. 
We must stress that our model is not Galilean invariant due to the form 
of velocity correlator ([3]) which reflects the existence of sweepping effect in 
developed turbulence [17]. As a consequence the renormahzation constants 
at ip'^dfip and ip'^{vd)ip terms are different in this case. Moreover, to assure 
multiplicative renormahzation of the model we have to add a new term ip'^ipvv 
into the action (jlj) with new independent parameter (charge) U2 which will 
compensate divergencies produced by the correlation function (ip^ipw) . The 
resulting renormalized action can be written in the following form 

= ^\-Zidt-Z2{v.V) + Z3DV^ -ZiDr)tlj 

_ ^ ^g^^t^v2 _ ivD-V. (7) 
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We note that due to the dimensional reasons the charge U2 must appear com- 
bination U2/D, otherwise it would possess nonzero dimension with respect 
to the spatial and time variable respectively. On the other hand the renor- 
malized action can be obtained by the multiplicative renormahzation of the 
fields ip — )■ ipZ^, ip"^ — ip^Z^p] and v — )■ vZ^ and parameters: 

Do = DZd, gio = gif^'^^'Zg^, to = tZ^, uiq = Ui^^'Zu^, 

Ao = Xn^Zx, 920 = g2f^'^^Zg2, U20 = M2^«2- (8) 

The renormahzation constant of fields and parameters ([8]) are related with 
renormahzation constant of the action ([7]) through the following relations 



Zv 


— Z^ ^Z2, 


Zr^ 




Zx 


= z- 


1/2 7-1 7 
Z5, 


Zd 


= Zi ^z^, 


Zux 


= ZiZ^ ^, 


Zr 


= ^3 




^91 


ryh 7 — 2 '7~'i 
— ^1^2 ^3 ' 


^92 


rv^2 ry2 


Zu2 


= Z2 


'^Z-j,Zi^. 



(9) 

The standard RG perturbative approach is based on the Feynman diagram- 
matic technique (TTJIIS]- Explicit calculation in the one-loop approximation 
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leads to following results for renormalization constants 



85' 85' 165 8^(1 



uje 



45 25 165 8mi(1 + Ui)e 

The basic RG differential equation for the renormalized Green function Gr 
is given by the equation {Drg + N^-^^ + A^^t7^/.t}Gij(e, /i, . . . ) = 0, where 
e is the full set of renormalized counterparts of the bare parameters Cq = 
{-Dq! To; ""lO! ^20) fi'iO) to} ^"^^ ■■■ dcuotcs othcr parameters, such as spatial 
or time variables. The RG operator Df^c can be written in the form 

Drg = /i<9^ + + /3«2^«2 + (^gidgi + f^g^dg^ - IdDd - 7r^r, (11) 

where = xd^ for any variable x, '-fa = In Za is an anomalous dimension 
and denotes the differential operator taken at fixed values of the bare pa- 
rameters. The P function is defined as Pg = D^g, g G {gi, g2, ui, U2} and from 
the relations (|8]) the following expressions can be obtained in straightforward 
manner 

^gi = ^i(-2e -V + 37d), Pg, = 92{-2S - 7gJ, 
Pui = Ui{-V + Id), Pu2 = -U2lu2- 

From the explicit results (fTOj) and relations ([9]) the needed anomalous dimen- 
sions are easily calculated 



35-1 , g2 ^9i ^9 



2 



4mi(1 + mi) 8' 2ui{l + ui) 2' 

92 , 3^1(1 + ^2) , . 

= + Au,il + u,y ^^^^ 



4 Fixed Points and Scaling Regimes 

According to the renomalization group theory, the infrared (IR) asymptotic 
behavior is governed by IR attractive fixed points (FPs). The fixed points 
9* — {5'*) fi'25 ^2} can be found from the requirement that all (3 functions 
simultaneously vanish f3g^{g*) = f3g2{g*) = (3^(9*) = f3u2{9*) = 0. The type 
of FP is determined by the eigenvalues of the matrix Q = {Q^ = df^i/dgj}, 
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where /3j is the full set of /3 functions ( TT^ and gj is the full set of charges 
{5'i7 5'2; Wi, M2}. For the IR attractive FP the eigenvalues of the matrix Q are 
strictly positive quantities. From this condition the region of stability for 
the given FP can be determined. It is well known that in IR asymptote the 
Green functions exhibit scaling behavior. The critical dimension Ap oi the 
IR quantity is determined by the following relation Ap = dp + A^dp + 7^,, 
where d^'^ are canonical dimension of quantity F from Tab. [T], 7^, is value 
of the anomalous dimension at the given fixed point and A^^ = 2 — 7I). 

Now we introduce the quantities [1], that could be used for illumination 
of the fixed points' structure of the model. The first quantity is the number 
of active particles N{t), which are spreading from origin x = at time t = 
and is given by the response function as follows 

N{t)= / dx(7/'(t,x)V^^(0,0)) oct^^ =t^^^^, (13) 



where 7I) = jD^gl, g2,u*). The second quantity is the survival probability 
P(t) of an active cluster and is defined by the following expression 
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P{t) = - lim (e-^^^O'^'V^l-^, 0)) oc = t . (14) 



A— >oo 



The radius of gyration R{t) of the active particles is defined as follows 

/ dx x2(^(t,x)v:'t(o,o)) 



R\t) = — ;y ' /;/ \ ' oct^ = (15) 

^ ^ /dx (^(t,x)V;t(0,O)) ^ ^ 

In order to investigate the FPs of this model, we use the exact relation for 
beta function I3g^/gi — 3(3uj^/ui = 2{rj — e). It follows that beta functions 
fig^ and I3ui can vanish simultaneously for finite value of their arguments 
only in the case e = rj, which should be studied separately. On the other 
hand, when e 7^ 77 it is necessary to consider limiting behavior, i.e. either 
Ml = or Ml = 00 and rescale gi so that in beta functions (3g-^, P^i anomalous 
dimension 7/5 acquires finite value. Let us consider case, for which limit of 
correlator velocity field is known as a 'frozen' velocity field and exhibits four 
FPs Tab. [21 In this regime the new variable g" = gi/ui was introduced and 
its beta function has the form Pg/^ = g'l{—2e + 27/)). The FP 1 correspond 
to the Gaussian (free) fixed point. In this regime the radius of gyration R{t) 
imitates the ordinary random walk (-R(t) oc t^^"^) behavior. For the FP 2 
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Table 2: Stable FPs of the "frozen" velocity field with corresponding critical 
exponents. 



FP 


FP 1 


FP 2 


FP 3 


FP 4a 


FP 4b 


gl" 








4e 
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4(6£"<5) 
15 


4(6e-<5) 
15 


g*2 





4<5 
3 





8{5-e) 
5 


8{5-e) 
5 





















not fixed 





-1 





25-le 
6e-<5 


Region 


e < 


6e<5 


e > 


< e 


< e 


of stability 


6 <0 


6>0 


e > 6 


|e < 5 < 6e 


e < 5 < |e 




T] <0 


6r] < 6 


Tj < e 


?7 < e 


1] < e 




2 


2- ^ 

^ 6 


2 - e 


2-e 


2 - e 




2 


7<5-12 
12-5 


5-2 
2-t 


6(5-e-10 
5(2-6) 


6<5-e-10 
5(2-e) 


{d - 2A^)/A^ 





25 

12-5 





2 (5-e 
5 2-e 


2 5-e 
5 2-e 



the correlator of the velocity field is irrelevant and model behaves like the 
"pure" DP class with corresponding critical exponents [1]. For the FP 3, the 
percolation nonlinearity of the DP action ([1]) is irrelevant. Finally, the FP 
4a and FP 4b corresponds to the nontrivial IR scaling regime, in which the 
nonlinearity of DP model ([T]) and velocity field are both important. Moreover, 
for the FP 4a the interaction i/j^ipv"^ is irrelevant unlike for the FP 4b. In the 
interval e G (0, 2) and for 5 =const. the exponent of survival probability P{t) 
is less than —1 for the FP 3 and greater then —1 for FP 4. It follows that 
in regime FP 4 the survival probability decrease more slowly than for FP 3. 
In the regime FP 3 the velocity fluctuations leads to the decay of an active 
cluster, but in regime FP 4 leads to the oposite effect. In the Fig. [1] the 
regions of stability for fixed point in the plane e — 6 for t] < are depicted. 
The boundaries of the regions are represented by the thick lines. 

The second limit case is characterized by white-in-time nature of velocity 
field. This regime is called the rapid change model [H] and exhibits four 
FPs. For this case it is advantageous to introduce new variables g[ = gi/u\ 
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Figure 1: Region of stability for fixed points of frozen velocity field. 



and w = 1/ui, for which corresponding beta functions have the following 
form I3g>^ = g[{ri — 2e + 'Jd), Pw = wij] — 7/)). The FPs and their regions of 
stability were published in the work Antonov et al. 



Table 3: Fixed points for the non-trivial case. 



FP 


FP 5 


FP 6a 


FP 6b 




4e 


4(6e-5) 


4(6e-<5) 


«*(!+«•) 


3 


15 


15 


9l 





8{5-e) 
5 


8(5-e) 
5 


u\ 


-1 





2<5-7e 
6e-<5 


Region 


e > 


e > 


e > 


of stability 


?7 = e > 5 


r/ = e<f 


r/ = e>f 




2-e 


2 - e 


2-e 




<5-2 
2-e 


65-e-lO 
5{2-e) 


6(5-e-10 
5{2-e) 


((i-2A^)/A^ 





2 <5-e 
5 2-e 


2 <5-e 
5 2-e 



The third, most non-trivial case, differs from the previous ones, the beta 
functions and ^5^^ are proportional to each other and thus the FPs are 
degenerated. However, the FPs can be divided according the value of 
into the two groups (see Tab. [3]). For the FP 5 the nonlinearity of the DP 
action ([T]) is irrelevant. For the FP 6a and FP 66 the interaction part of DP 
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action ([T]) and correlator velocity field coexist. In the interval e G (0,2) the 
exponents of quantities {N(t), P(t), R{t)) behave in the same way as FP 3, 
FP 4a and FP 46, respectively. 

5 Conclusion 

In this work we have investigated the influence of the finite correlated velocity 
field (Antonov-Kraichnan model) on the progress of directed bond percola- 
tion of process. After the construction of the field-theoretic model, we have 
performed the calculations of the renormalization constants of the triple ex- 
pansion in e, 2(5 = 4 — (i and t] in the leading order. Fixed points, their 
regions of stability and behavior of the number of active particles, survival 
probability and radius of gyration were calculated. 
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